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Abstract 

In this series of lectures we present the universal method based on 
the use of the functional integrals to derive the bound state equations 
for different two-body systems in elementary particle physics as well 
as in condensed matter theory: the relativistic Nambu-Jona-Lasinio 
equation and the relativistic Bethe-Salpeter equation in elementary 
particle physics, the Bethe-Salpeter equation for the bipolaritons. 
The bound state equation for the three-body systems is also estab- 
lished. 

1 Introduction 

The theoretical study of the bound states of different systems of fundamental particles, 
their energy (mass) spectra and the corresponding many-body quantum-mechanical 
wave functions or quantum field-theoretical state vectors was always one among the 
main problems of the quantum physics. The simplest examples of the two-body quan- 
tum systems are the hydrogen atom and the positronium. In the non-relativistic 
approximation for the proton the energy spectrum and the wave functions of the hy- 
drogen atom are determined by the Dirac equation with the Coulomb potential of the 
attractive interaction between the electron and the proton. For the formation of the 
bound states of two particles (including particle-antiparticle pairs) due to the boson 
exchange in the relativistic quantum field theory, for example the positronium in QED, 
we have the Bethe-Salpeter equation' 1 !. This relativistic integral equation was widely 
applied to study the bound states of different two-body systems: the mesons as the 
bound states of the quark-antiquark pairs in QCD' 2 ~ 4 1 and the diquarks as the bound 
states of the quark-quark pairs in dense QCD' 5 ' 6 '. The direct four-fermion coupling 
may be also the physical origin of the formation of the composite particles as the 
bound states of some fermion-antifermion pairs or two-fermion systems. In this case 
we have the Nambu-Jona-Lasinio equation' 7 !. It was widely applied to the study of 
the formation and the physical properties of the mesons as the bound states of the 
quark-antiquark pairs' 8 " 15 ! . Recently this equation was used for the investigation 
of the instanton induced quark-quark pairing in dense QCD' 16 ~ 21 1. In the study of 
the formation of the bound states of two-body systems of quasiparticlcs in condensed 
matters one can use the Schrodinger equation if the effective mass approximation is 
valid. However, there are quasiparticles with the complicated energy spectra so that 
we cannot establish the Schrodinger equation. In this case we must derive and then 
apply either the Bethe-Salpeter or the Nambu-Jona-Lasinio equation. For example, 
due to the exciton-exciton interaction there may exist the bound states of the sys- 
tems of two polaritons called the bipolaritons. The Bethe-Salpeter equation for the 
bipolaritons was established and studied by many authors' 22 ' 23 !. 

The relativistic Bethe-Salpeter equation for the systems of two elementary particles 
as well as the nonrelativistic Bethe-Salpeter equation for the systems of two quasipar- 
ticles with the complicated structure of the energy spectra in condensed matters were 
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derived within the framework of the perturbation theory by summing up the matrix 
elements of the infinite series of Feynman diagrams in the ladder approximation. The 
Nambu-Jona-Lasinio equation was established in the pioneering work of Nambu and 
Jona-Lasinio also by means of the perturbation theory! 7 ! . Subsequently the elegant 
derivation of this famous equation by means of the functional integral technique was 
given by Eguchi, Sugawara' 8 ' 9 ! and Kikkawa' 10 !. Within the functional integral for- 
malism it is not necessary to sum up the infinite series of the matrix elements of the 
Feynman diagrams in the ladder approximation, and the derivation of the equations is 
simplified. In the attempt to present the functional integral technique as the universal 
method for the study of many different problems in the quantum theory! 24-29 ! in this 
scries of lectures we review the derivation of both types of bound state equations in 
the elementary particle physics as well as in the condensed matter theory by means of 
the functional integral technique. The derivation of the bound state equations for the 
systems of two quasiparticles in condensed matters by means of the functional integral 
method is presented first time in these lectures. As we shall see, it is straightforward 
to establish the bound state equations for the three-body systems' 30 ' in the framework 
of the functional integral formalism. 

2 Relativistic Nambu-Jona-Lasinio Equation 

We start to study the equations for the bound states of many-body systems by con- 
sidering the simplest examples : the formation of the mesons and the diquarks due to 
the direct 4-fermion coupling of the quark field with the interaction Lagrangian 

L int (*) - ^W^W^VW^W, (1) 

tjBD _ ttBD _ tjDB _ ttDB 
U AC — ~ U CA — ~ u AC — U CA i 

where tjj A (x) denote the quark field, 

A = {aai) , B = (fibj) , C = (jck) , D = (Sdl) 

are the sets consisting of the Dirac spinor indices a, (3,j,S = 1, 2, 3, 4, the indices of the 
color symmetry a, b,c,d — 1, 2, ...N c and those of the flavor one k,l — 1, 2, ...Nf. 
For the convenience in the study of the diquark formation we write the interaction 
Lagrangian (1) also in another form 

Lint (X) = (x) ~f (x) V^A^D (*) ^1>B (*) , (2) 

with the new coupling constants 

vEZ - KB- (3) 

The mathematical tool in the functional integral formalism is the functional integral 
of the field system 

Z = J [Di>] [L^]exp{iS[V,?]}, (4) 
where S [ip, ^T] is the total action of the system 

S[iP,^]=S ty,J\+S int [iP,J\, (5) 
So [tp, ^] being the action of the free quark field with some bare mass m 



B ri 
+ TOO, 

A 



ip B (x) , (6) 
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and Sint [V>, W\ being the contribution of the interaction Lagrangian to the total action 

S^t [V,?] = f dtxL^ix). (7) 



For the study of the mesons as the bound states of the quark-antiquark pairs we use 

the intei 
and set 



the interaction Lagrangian in the form (1), introduce the composite meson field $^ (a;) 



Zt = J [D$]exp{-1 J d^ A B (x)U*°Z% (*)}. 

Shifting the functional integration variables 

(z)-^ (x)+^ A (x)^ B (x), 
we establish the Hubbard-Stratonovich transformation 



(8) 



exp | l - J d^ A (x) *P B (x) U B A S^ C (x) 4> D (x) | = (9) 
± J [D$]exp|-i J d^i (x)USB^ (*)}exp{-i J d^ A (x)^ B (x) A B 



with 



A%(x) = U%£$ D {(x) (10) 



Using this formula, we transform the expression (4) of the functional integral of the 
quark field into that of the composite meson one 

z = %\ m]^p{is eS m} (ii) 

with the effective action 

S cS m = ~\j (x) U%8* C d (*) + W [A] , (12) 

where 

Z = J m\ [D$\ exp {iS [V,?] } (13) 
and the functional W [A] is determined by the formula 

exp {iW [A] } - i- /" [£ty] [ity] exp {zSq [V, } 

exp |-*y rf4 ^ 0*0 1>b {x) Af (a:) | . (14) 
We can express W [A] in the form of a functional power series in the field A^ (a;) 

oo 

W [A] = ^2 W {n) [A] , (15) 



n=l 



where [A] is a homogeneous functional of the n-th order. Calculations give, for 

examples, 

W^[A] = -iJd i xS^(0)^i(x), (16) 
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[A] = l - J d 4 x J d 4 yS^ [x y) Ag (y) Sg (y x) A£ (x) , (17) 



[A] = -i y d 4 x I d 4 y | d 4 zSf (a; - y) Ag (y) S£ (y z) Ag (z) S| (z x) 



(18) 



etc, where S B (x — y) is the two-point Green function of the free quark field 



Sa (x-y) = 



Z 



J [D$\ [D$\ ip a (x) ^ (y) exp { t S [V,?] } • (19) 



It satisfies the equation 



d ) + m5% 



S c B {x-y)=5 c A 5(x-y). 



Denote S B (p) its Fourier transform 

S a^) = 7^TS / d 4 pe^S B (p). 



We have 



1 



From the variational principle 



8<S>% (x) 



ip + m 



= 



we derive the field equation 



A c (x)-C/ CA ^ 



(20) 



(21) 



(22) 



(23) 



(24) 



For some definite types of composite meson fields the first order functional [A] either 
vanishes or does not give the contribution. Then in the second order with respect to 
A^ (a;) we can replace W [A] by [A] and obtain the approximate field equation 



A c {x)-U CA SA B [x) 



(25) 



Substituting the expression (17) of [A] into the r.h.s. of the relation (25), we 
obtain immediately the Nambu-Jona-Lasinio equation 



AS (x) = zUSi 3 J d'ySE (x y) A F E (y) S$ (y - x) . 



(26) 



For establishing the field equation in the general case with non-vanishing [A] we 
must take into account the contributions of the functional [A] of all orders by 

summing up the infinite series (15). Denote 

S B A {x,y) = S B (x-y)- Jd 4 x 1 S^(x-x 1 )Ail(x 1 )Sl(x 1 -y) 

+ J d 4 x Y J d 4 x 2 Sl (x - Xl ) A* ( Xl ) Sfl ( Xl - x 2 ) Ail O^) Si (x 2 - y) 

- J d 4 Xl J d 4 x 2 J d 4 x 3 S Bl {x - Xl ) A B \ ( Xl ) S B l ( Xl - x 2 ) 

(x 2 ) S B A l (x 2 - x 3 ) A^ 0*0 Si (x 3 -y) + (27) 
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the two-point Green function of the quark field interacting with the composite meson 
field (x) . It is determined by the Schwinger-Dyson equation 

Si (x, y ) = S %{x-y)- J d'zS c A (x - z) Ag (z) Sg (z, y) . (28) 

In terms of this new two-point Green function we can rewrite the field equation (24) 
in the form 

AS(x) +i Ug*Si(x,x)=0. (29) 

In general there exists some non- vanishing constant solution of the nonlinear equation 
(29) 

(<E> )^ = const, 

(A )* = U^^°f D = const. (30) 

With the constant meson field the expression in the r.h.s. of the formula (27) depends 
only on the difference x — y of the coordinates x and y 

Si(x,y) = Si(x-y), (31) 
and this constant meson field is determined by the equation 

(A°)2+iC#Ts2(0) = 0. (32) 

Consider the quantum fluctuations of the meson field around the background constant 
one (30) satisfying the equation (32) and set 

Ai(x) = (A°) A B +^(x), (33) 

with 

it (x) = U£g<pS (x) . (34) 

Up to the second order with respect to the quantum fluctuations the effective action 
of the composite meson field equals 

s cS [*] « s cS [*°] + \ J d-x J dSA (x) s ^to S *c iy) \~*<fZ (y) ■ ( 35 ) 

Using the expression (12) of S e f[ [<&] , we have 

+ l/^/^M , A /^i fe , L^£(»). (36) 

By summing up the infinite series 

S 2 W^ [A] 



^SAi (x)5A% (y) 
we obtain 

S 2 W[A) , c 



a=aO 



8A* (x) 5A% (j/) Ia=a 



=iS%(x-y)S£(y-x). (37) 
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Then from the formula (36) for the effective action and the variational principle it 
follows the Nambu-Jona-Lasinio equation containing the two-point Green function of 
the quarks interacting with the constant background meson field 

d (*) = iU%8 J d 4 yS% (x y) f E (y) S C F (y x) . (38) 

In the special case when the background meson field vanishes or its contribution is 
neglected then the equation (38) reduces to the equation (26). 

Now we study the formation of the bound states of a pair of two quarks. For this 
purpose we use the interaction Lagrangian in the form (2), introduce the composite 

bosonic diquark field $ab (x) as well as its conjugate $ (x) and set 

Z** = j [£>*] [DW] cxp j-1 j d^ AC (x) V$?* DB (or) J . (39) 
Shifting the functional integration variables 



$db(x) -> $ DB (x)+4> D (x)ip B (x), 

¥ AC (x) -» ¥ AC ( x ) + i> A {x)^ C (x) , (40) 

we establish the Hubbard-Stratonovich transformation 



\ J d 4 x^ A (x) (x) V^ D (x) ^ B (x) 
1= J [£>$] [£>$] exp |-1 1 (*) V^f <& DB (*)} (41) 

| - \ j d A x$ A {x) ~f (x) A c a {x) + A BD (x) yj D (x) 4> B (x) J , 



Z 

. exp 



where 



A CA (x) = Vi^ DB {x), 

A BD {x) = $ AC (x)V c B £. (42) 



Using formula (41), we can transform the expression (4) of the functional integral of 
the quark field into that of the diquark system 



Z = -|= / [D*\ exp {iS eS [$, $] } (43) 
with the effective action 

S cS [$,¥] =~J d A x<§ AC (x) V§2* DB (x) + W [A, A] , (44) 
W [A, A] being determined by the formula 

exp {iW [A, A] } = / [Ity] exp {zS ty, ^] } (45) 



CX P ^ ~ 2 ^ ^ ^ ^ ^ C ^ AcA ^ + ^ ^ D ( X ^ B ^ 



and having the form 
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W [A, A] = ^2 [A, A] , (46) 



71=1 



where W^ 2 ^ [A, A] is a homogeneous functional of the n-th order with respect to each 
type of bosonic fields Aq A (x) and A (x) . Calculations give 

[A, A] = ~J d A x, J d 4 x 2 A AlBl ( Xl ) S% {xi - x 2 ) S*\ (an - x 2 ) Ab 2 a 2 (x 2 ) 

(47) 

or in another form 

[A, A] =~J d 4 Xl J d 4 x 2 A AlBl ( Xl ) Sil ( Xl - x 2 ) A B2A2 (x 2 ) S^ 2 (x 2 - Xl ) 

(48) 

convenient for the comparison with the higher order functionals 

[A, A] = l - J d 4 Xl ...J d 4 Xi A AlBl ( Xl ) Sil (x t - x 2 ) A B2 a 2 (x 2 ) Slf 2 (x 2 - x 3 ) 
A A3B3 (x 3 ) S^l (x 3 - x A ) A BiAi (x 4 ) S T A Ai (x 4 - Xl ) , (49) 

[A, A] = - l - J d 4 Xl ...J d 4 x 6 A AlBl (x 1 )S^ i (x 1 -x 2 )A B2A2 (x 2 )Sl^(x 2 -x 3 ) 

A MB " (x 5 ) S§1 (x B - x 6 ) A BeAe (x 6 ) Slf e (a* - xi) . (50) 

etc, where we used the notation 

SJ B (y -x) =Si(x-y). (51) 
From the variational principle 

5S cS [$,¥] 



OK™ ^ = ° (52) 
8$ (x) 



we derive the field equation 



fato.VgnfZ&S.. (53) 
1 SA (x) 

If we neglect the contributions of the high order functionals W^ 2n ^ [A, A] with n > 1 
and use the approximation 

W [A, A] « W {2) [A, A] , 

then the field equation (53) is the Nambu-Jona-Lasinio equation for the composite 
diquark field 

A AB (x) = -iV»§ J d 4 ySg (x - y) S£ (x - y) A EF (y) . (54) 

In order to take into account the contributions of the high order functionals W^ 2 ^ [A, A] 
with n > 1 we search the non- vanishing constant solution of the field equation (53) 

AC 



(® ) DB = const, (<I> J = const, (55) 
(A°) CA = Vi2{&) DB = const, (A°) BZ5 =(¥°)' 4C ^=const 
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and consider the quantum fluctuations around this background constant field 



<P DB (x) = (*°) DB +<p DB (x), $ AC (x) = ($°) AC + ^(x), 
A CA {x) = (A°) CA +t CA (x), A BD (x) = (^) BD +f D (x),(56) 
with 

tcA {x) = VSaVdb {x) , f° (x) = W AC (x) V^f. (57) 
Up to the first order with respect to each type of bosonic fluctuations <p DB (x) , £ >DB (x) 
and^ c (x),C (x) we have the effective action 



S cS [$, *] « S cB [*°. * J = 5 / d4xJ P AC ( x ) Obb (*) (58) 

+ fd 4 x fd 4 yl AC (x) Ac^Jf'^ , H a=ao,a=a°^b (*) • 
<5A (x) (y) 

By summing up the infinite series 

8 2 W [A, A] 



we obtain 



<y 2 w [A, A] 



n= i (5 A (x) 5A DB (y) 



^ac, rrr — ia=ao,a=a- = 4 s * (* - y) sg (z - y) , (59) 

dA (x) SA DB (y) z 
where (x — y) is determined by the Schwinger-Dyson equation 

(x-y) = Si (x y)-j d 4 z J d A uS c A (x - z) (A a ) CD S T E D (z u) (a )^ Sf (u-y). 

(60) 

Then from the expression (58) of the effective action and the variational principle it 
follows the Nambu-Jona-Lasinio equation 

Ub (x) = -iVRE I d 4 yS% (x - y) S£ (x - y) £ EF (y) (61) 

containing the two-point Green function (a; — y) of the quark field interacting with 
the background bosonic constant field (55). In the special case of the vanishing back- 

/ q\ AC / q\ AC 

ground bosonic field ($°) DB = ( $ J = (A°) DB = ( A J = the equation (61) 
coincides with the equation (54). 

In order to solve equation (60) we work in the momentum space. Denote S A (p) 
and (p) the Fourier transforms of the two-point Green functions S A (x) and (x) , 

Si(x) = ^4 / d*pe^S B A (p), (62) 
px — px — pox = px Et, 
and introduce the matrices S (p) , S T (—p) , S (p) , A and A° with the elements S A (p) , 

~ ~. / \ AB 

S A (—p) , S (p) (A ) AB and A . Then we can rewrite the Schwinger-Dyson equa- 
tion (60) in the matrix form 



8 



1 



1 



A°S T (~p)A^. 



(63) 



S(p) S(p) 

For the constant background composite field (55) the field equation (53) has the explicit 
form 

(A°) AB = -iVft {j d±ySB( X -y)(^) EF Sl F {y~ X ) 

-Jd 4 yJ d 4 x, J d 4 Vl S F (x y) (A ) EF S T D F (y Xl ) 

. (A°) DlCl S% ( Xl - Vl ) (A°) EiFi ( Vl x) (64) 

+ J d A yJ d 4 Xl J Ai / d A x 2 J d%S F (x - y) (A°) EF S T D F (y - x x ) 

. (A )" 1 " 1 S% ( Xl - m ) (A°) ElFl Si? ( Vl - x 2 ) (a°)^ 
.S*l{x 2 -y 2 )(A*) E2F2 S T D F >{y 2 -x)- }. 

In term of the Fourier transforms of the two-point Green functions this equation be- 
comes 



(A°) 



AB = -iV^^s J d i p{sE{p){A°) EF S F D {-p) 



(2nf 

-SB (p) (A°) £F S F Dl (-p) (A*) DlCl S F l (p) (A°) EiFi S* (-p) 
+S F (p) (A°) EF S F Dl (-„) (A°) DlCl S F l (p) (A°) EiFi S F l (-„) 
.(a ) D2C2 SE : (p)(A°) E2F J f H- p )- }. (65) 



or in the matrix form 



1 ,AB AB (2n) 4 J *|_ l l + A^S{p)A»ST{-p)\ CD 

(66) 



3 Relativistic Bethe-Salpeter Equation 

Now we generalize the reasonings presented in the preceding Section to establish the 
bound state equations for the two-body systems of relativistic Dirac fermions with some 
effective (non-local, in general) 4— fcrmion interaction in the quantum field theory For 
the definitcncss wc consider again the quark field with some effective direct 4— quark 
(non-local, in general) interaction induced by the dynamical mechanisms in QCD. 

Introduce the Fourier transforms ip A (p) and i/j (p) of the quark field ip A (x) and its 
conjugate tp (x) , 



i>A 0) = 



ip A (x) 



i-2 J d'pe^A (P) , 



(27T) 
1 

(2^ 



d 4 pe- ipx ip (p). 



(67) 



Then we have 
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and 



S [l>,tf\=-J d 4 P ?P (p) [i (p) B A + m6%] tpg (p) (68) 

J- J [Dip] [Dip] i> A (p)f (q) cxp {iS [4>,tf\} = 6(p-q) S* (p) (69) 

The contribution Si nt [ip, ^TJ of the 4— quark interaction to the effective action repro- 
ducing the 4— point Green function of the quark field will be represented either in the 
form 

Sint = \j d i p 1 d i q 1 J d A p 2 d A q 2 5^ (Pi + qi- P2 - <h) (70) 

~M ~ A ~A 2 

4> (Pl)V'Bi (-<H) U a1b2 (PU-qi; -Q2,P2)lp (-q2)lp B2 (P2) 

convenient for the study of the mesons as the bound states of the quark-antiquark 
pairs, or in another form 

Sint - \ j d A p 1 d A q 1 J d A p 2 d A q 2 S (A) (p! + qi -p 2 - q2) (71) 

V> (Pi)V> (<?i) V^^ 2 (qi,Pi;P2,q2)4>B 2 0?2)V% O2) 

convenient for the study of the bound states of two quarks - the diquarks. The covariant 
quantities U A * A ^ (pi, qi;P2, 92) and V^ 2 ^ 2 (<7i,Pi, ;P2, 92) may be called the relativistic 
effective potentials. They arc related each to other 



^rir (01 , pi ; w , <&) = ttf (pi , p 2 ; qi , g 2 ) (72) 

and must be antisymmetric under the simultaneous interchanges of the upper or lower 
spinor indices and the corresponding 4— momenta of the quark field 

u mbI (pi>-9i;-?2,P2) = - U^IaI (-52,-gi;pi,92) = 
= -^b, 1 (pi>P2; -92, -91) = tf^ 1 (-<Z2,p 2 ;pi, -91) , 

V B 4 1 2 i ? 1 2 (9i,Pi;P2,g 2 ) - -V£ a * a (pi,gi;p2,<fe) = (73) 
= -^aI 2 Pi! 92, P2) - ^ 2 (pi, gi; ga.pa) • 

In order to study the mesons wc start from the expression (70) of Si nt 
introduce the composite bi-local meson field (x, y) , whose Fourier transform is 
denoted (p, q) 

$b {x, V) = ^4 / d A pd A qe-^-^l A B (p, q) , (74) 

and set 

Zt = J [£>$] cxp |-1 y <i 4 pi Al y d 4 P 2d 4 q2S (4) (pi + gi - p 2 - 92) 

fa,-<h)U%£ (Pi,-9i;-92,P2)$f 2 (-92, P2)} (75) 
Shifting the functional integration variables 
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<f b (p, -9) - (P, -9) + / (p) V' s (-9) » (76) 
we establish the Hubbard-Stratonovich transformation 

CXP { \ J d4pid4qi J rf4 P2rf 4 92<5 (4) (pi + qi - p 2 - 92) 

^ (Pl)^B! (-9l)^AiB2 (Pl,-'7i;-'72,P2)V' S2 (-92) ?/>A 2 (P2)| = 
= ^* / ^ ° XP / dipidiqi j d4 P2^ 4 92^ (4) (Pi + 91 -P2 - 92) 

(pi,-Qi)U%£ (Pi,-?i;-92,p 2 )^ (-<Z2,P2)} (77) 
expj-i y d 4 pd 4 qi{> (p) V> 5 (-?) Af (p, -q>) 



where 



A* (p, -«) - I d 4 p'dWS {4) (p + q-p'- q') U B J, (p, -q; -q',p') (-q',p') . 

(78) 

Applying the transformation (77) to the expression in the r.h.s of the formula (4) 
with the effective interaction action 5j nt of the form (70) and performing the 

functional integration over the fermionic integration variables ip A (x) and ip A (x) or 



-A 



tp A (p) and ip (p) , we rewrite the functional integral (4) of the interacting quark field 
in the form of the functional integral (11) of the composite bi-local meson field with 
the effective action 

S eS [$] = -\\ Ai Ai / d i p 2 d i q 2 8 {4) (pj + qi - P2 - 92) (79) 
( Pl , - gi ) (pi, -gi; -g 2 ,p 2 ) $t (-92,P2) +W[A], 

where W [A] is a functional power series of the form (15) in the bi-local meson field. 
Calculations give 

[A] = -iJ d'pSi (p) A* (p, -p) , (80) 

[A] = 1 I dV I rf 4 p 2 ^ (pO A* (p 1; -p 2 ) (p 2 ) A* (p 2 , - Pl ) , 

(81) 

^ (3) [A] = - l - J d A Pl J d A P2 J d 4 P3 s£l (pO A* (p l5 -p 2 ) 

S£ (pa) Af« (p 2 , -p 3 ) S£ (P3) (P3.Pi) . (82) 



etc. 

From the variational principle 



<S$£ (p, -g) 
it follows the field equation 



5SM =0 (83) 
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A B (p, -q) = J d 4 p'd 4 q'5^ (p + q-p'- q') U B J, (p, -q; q',p> 



SW [A] 



^'(-^p') 
(84) 

If [A] vanishes or docs not give the contribution, then in the second order approx- 
imation the field equation (84) becomes a homogeneous linear integral equation for the 
bi-local meson field - the relativistic Bethe-Salpeter equation for the quark-antiquark 
bound states 

Af (p, -q) = i j d 4 p'd 4 q'S^ (p + q-p'- q') 

UIb> (p, ~Q\ -q',p') S% (p') Ag (p', -q>) S B ' (q') . (85) 

In the general case, when [A] is non- vanishing or the contributions of the high- 
order functionals [A] , n > 1, are not negligible, we start to search the solution 
of the field equation (84) in the class of the function of the special form 



(*°(p,-g)) B = 5^(p + q)<S>i(p), [A (p,-q)) A =S^(p + q)^(p), 
Af (p) = J d 4 p'U BB ' (p, -p; p', -p') {-pf) , (86) 

which are the Fourier transforms of the functions ($° (x — y))^ and (A (x — y)) A 
depending only on the difference x — y of the coordinates, 

(*° (x - y)) A B = ^4 / rfV 4 9^' (pl+? ^ (4) (p + q) H (P) 

= (2^ / di P e ~ V{X+V) ®B(p), (87) 
(A (x y)) B A = ^4 / dW- <(px+w) * (4) (P + ?) Af (p) 

-i-j [dtpe-W^Aiip), 
(2,1:) J 

and then to consider the quantum fluctuations of the bi-local meson field around the 
background one of the form (87) 

H(P,-q) = S^(p-q)^(p)+^(p,-q), 

A B (p,-q) = 5^(p-q)A B (p)+f A ( Pl -q), (88) 

Ta (P, -q) = J d 4 p'd 4 q'S {4) (p + q-p'- q') U B g, (p, -q; -q',p')f A , (-q',p') . 

(89) 

For the background meson field (87) the equation (84) becomes 

A B (p) = -ij d 4 qU B S (p, -p; -q, q) §g (q) , (90) 

where (p) is the Fourier transform of the two-point Green function of the quarks 
interacting with the meson field (86) : 
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Sf(p) = S^(p)-S^(p)A^ i (p)Si i (p) 

+S£ (p) ~A B A \ (p) Sil (p) All (P) Si 2 (p) - (91) 

It satisfies the Schwinger-Dyson equation 

§f (p) - (p) - ^ (p) Ag (p) Sg (p) . (92) 

In the second order with respect to the quantum fluctuations the effective action (79) 
equals 

5* cff m = S* cff [*°] - 1 y d*pi Ai y d i p 2 d 4 q 2 S^ (pi +p 2 - «i - ga) 
<Pb\ (Pu -Qi) Ua\bI (Pu -3i ; -<?2,P2) ^ (-92, p 2 ) + 
+1 y d 4 pi Ai y A2 A2 • (93) 

~A! 5 2 I^ [A] ~A 2 

5A Al (Pl.-9l)^ A A2 (P2.-92) ^ ^ 

By summing up the infinite series 

S 2 W^ [A] 



S^fr.-aOfA* (P2,- 9 2) IA -( A0 ^ 
we obtain 



r7 _Bi , ^^J-L , =(*.)- - ^ W (Pi + 92) 5 (4) (P2 + 9i) S* ( Pl ) S* (Pa 

SA A \ (p!,-gi)5A^ (p 2 ,-g 2 ) 1 ^ 

. (94) . 

Then from the formula (93) for the effective action and the variational principle it fol- 
lows the relativistic Bcthe-Salpeter equation containing the two-point Green function 
of the quarks interacting with the background meson field 



U (p, -q) = J d 4 p'd 4 q'S {4) (p + q-p'- q') 



U^l (p, -q; -q',p') S% (pO (p'> Sg' (</) • (95) 

In the special case when the background meson field vanishes or its contribution is 

neglected then the equation (95) reduces to the equation (85 ) 

In order to study the formation of the diquarks we start from the expression (71) 

of Si n t introduce the composite bi-local bi-spinor bosonic field &ab (x,y) and 

—ab ~ — AB 

its conjugate <I> (x, y) as well as their Fourier transforms &ab (p, q) and $ (p, q) 

*AB (x, y) = / AA^+^ab (p, q) , 

4R 1 f , ^_AB 

* (*,y) = -j rfVV" (pi+w) $ (p,«), (96) 

(27TJ J 

and set 

^ ,¥ = y cx P |-| y ^pi y Ai y a 2 y (pi ^ 



$ (Pl,9l) Vslll (9l,Pi;P2,«!) $B 2 A 2 (<ft,P 2 ) f> (97) 
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Shifting the functional integration variables 



®AB(p,q) ->■ $ab (p, q) + -4>a (p) (Q) i 
$ (p,g) - $ (p,?) + V (p)V> («), (98) 
we establish the Hubbard-Stratonovich transformation 



exp 1 1 y d 4 pi y d 4 gi y d 4 p 2 y d 4 <72<5 4 (pi + q\ - Pi - qi) 

i> (pi) V> (qi) v b ?m («i'Pi;P2, 92) V>s 2 (92) Vu 2 (P2) 

= -i= y [z?$] exp j-i y d 4 P1 y Ai y ^ y a 2 <* (4) ( P1 + qi - P2 - <») 

4 AlBl (Pi.gi) («i,Pi;P2,<B)*B a x a (?2,p 2 )J (99) 

(if f — A — B ~ ^_AB „ „ 

' exp 1 2 / rf4p / ^ ^ ^ AsA + A (P' q ^B (q) (p) 



&BA(q, P ) = J d A p' f d 4 q'6^ (p> + q> -p-q)V B A A B ' (q,p;p>,q')<l> B ,A'(q',p'), 

^_AB f f ^_A'B' 

A (p, q) = / dV / d\'5^ (p' + q'-p-q)^ (p' , q') V£? A , (q',p';p, q) . (100) 



The bi-local bi-spinor fields <&ab (p, q) and $ (p, q) are antisymmetric under the 
simultaneous permutations of the coordinates and the bi-spinor indices 



*ba {q, p) = -<$>ab (p, q) , Aba {q, p) = - A AB (p, q) , 

2^_BA ^_AB ^_BA ^_AB 

$ (q,p) = -$ (p,q), A (q,p) = -A (p,q). (101) 

Substituting the expression in the r.h.s. of the relation (99) into the functional integral 
(4) and performing the functional integration over the fermionic integration variables 

—A~— A 

ip A (x) and ip (x) or ip A (p) and ip (p) , we rewrite the functional integral (4) of the 
interacting quark field in the form of the functional integral (43) of the composite 
bi-local bi-spinor bosonic diquark field with the effective action 

S cS [$,¥] = ~\j d A Pi J Ai J d A P 2 J d%S 4 (pi + qi - P2 - q 2 ) (102) 

^ AlBl (Pi,9i)< 2 ^ 2 (qi,Pi;P2,q 2 )$B 2 A 2 (qi,Pi) + W [A, A] , 

where W [A, A] is a functional power series of the form (46) in the fields A A b {x, y) 

—ab ~ — AB 

and A (a;, y) whose Fourier transforms are Aab (p, q) and A (p, q) resp., 



Aab (x, y) = — L f d 4 pd 4 qe i(px+gv) A A B (p, q) , 
(2tt) J 

A (a;, y) = — I / d^qe^+^A (p,q). (103) 
(2tt) J 



Calculations give 
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[A, A] =- % -J d 4 p J d\k AB (p, q) SE (?) S C A (p) A DC (q,p) (104) 

[A, A] = 1 1 d 4 Pl J d% J d 4 p 2 J d%^ AlBl ( Pl , qi )S% (qi)A DlCl (qi,pi) 
S% (p 2 )t A2B2 {p 2 ,q 2 ) S% (q 2 ) A D2 c 2 (q 2 ,p 2 ) S% ( Pl ) , (105) 

[A, A] = l - J d 4 Pl J d 4 Ql J ...J d 4 p 3 J d 4 q^K lBl (pi,qi)S% (qi)A DlCl (?i,Pi) 
S% { P2 )1A A2B2 ( P2 ,q 2 ) SE 3 3 (<fe) Ad 3 C 3 (<a,P3) S% (px) , (106) 



etc. 

From the variational principle 



5S eS [$,¥1 

-AC = ( 10? ) 

(5$ (a;,y) 

it follows the field equation 

\a C a (q, P ) = I d 4 p' f d 4 q>5^ {p + q-p 1 - q 1 ) Vg% (p, g') ^j A ' A] ■ 
J J 6A {p>,q>) 

(108) 

In the second order approximation with respect to the composite diquark field we have 
the linear integral equation 

Aca (q, P) = -if d 4 p> J d 4 q'S^ (p + q-p' -q') 

v8a (p, q;p', q') s§ (?') s£ ( q ') a fe ( q ', P ') , (109) 

which is the relativistic Bethe-Salpctcr equation for the diquarks. In the general case, 
when the contributions of the high order functional W^ 2 ™) [A, A] , n > 1, cannot be 



neglected, we search the solution ($° (y — x)) DB ,l&(x — y)\ or (A (y — x)) CA , 
/— o \ Ac 

( A (x — y) J of the field equation in the class of the functions depending only on 

the difference x — y of the coordinates and study the quantum fluctuations around this 
background field: 



<S> DB (y,x) = ($° (y-x)) DB +ip DB {y,x), 
$ BD (x,y) = (*° (x-y)y D + tp BD (x,y) , 
A CA (y,x) = (A°(y-x)) CA + £ CA (y,x), 
A AC (x,y) = (A°(x-y)) AC +f C (x,y). (110) 
The Fourier transforms of the background field and its conjugate have special form 

q \ BD BD 

* (P,«)J =<5 (4) (p + <?)¥ (p), 

(_ q \ AO AO 

A (p,g)J = ^ 4 >(p + g)A (p), (111) 
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Denote (p DB (q,p) , ip (p, q) and £ CA (q,p) , £ (p, q) the Fourier transforms of the 
fields tp DB (y,x) , Tp BD (x,y) and ( Cyt £ ( x ,y), resp.. In the second order 

with respect to the quantum fluctuations the effective action equals 

s cS [$,$] « s cff [$°, $°] - 1 y d* P1 y Ai y a 2 y ^ + qi - P2 - q2} 



V 1 1 (Pi , 9i ) ^i 7 ! 2 (Pi > 9i ; P2 , 92 ) ^C 2 A 2 (92 , P2 ) 
+ / Al / Al / <**P2 / d% 



(112) 



-Aid 

€ (p,«) 



6 2 W 



^47cT- — " -\a ca = { a° )ca ,a ac =(a°) ac Zc 2 a 2 (92, pa) • 

oA (pi,gi)(5Ac 2 A 2 (92, P2) 



Summing up the infinite series 



A , A 



«=i<5A (p 1 ,q 1 )SA C2 A 2 {q2,P2) 



Ia ca = (AO) c ,A- 4C = (A°)^' 



we obtain 



where 



sS(p) 



A°,A U 



(5A (pi,gi)5Ac 2 A 2 (92, P2) 



l A C A = (AO) CA ,A =(A ) 



^ (4) (Pi - P2) <5 (4) (91 - 92) Sgj (px) , 



(113) 



5^ (p) - (p) A ClAl (p) s£ (-p) 3^' (p) £g 2 (p) 
+S£ (p) A ClAl (p) Sil (-p) 3^' (p) 5g3 (p) a C3A3 (p) 

<M-p)^ 4C4 (p)^ 4 (p)- 



(114) 



is the two-point Green function of the quark field in the presence of the background 
diquark one and satisfies the Schwinger-Dyson equation 

§S (P) = S C A (p) - S% (p) A ClAl (p) S£ (-p) A" A2C2 (p) Sg 2 (p) , (115) 
or in the matrix form 



S(p) 5(p) 



+ A(p)5 T (-p)A(p), 



(116) 



■AC 



where A (p) and A (p) are the matrices with the elements Ac a (p) and A (p) . 
From the expression (112) of the effective action and the relation (113) it follows the 
relativistic Bethe-Salpeter equation 



He a (9, p) = -i y A' y A'<5 (4) (p + ? - p' - ?o 

V£? (P, 9;p', q') Sf (p') (g') £ FB (9',P') 



(117) 
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containing the two-point Green function of the quark field interacting with the back- 
ground diquark one. If the later vanishes, then the equation (117) reduces to the 
equation (109). 

The relativistic effective potential U^g* C^ 1 ' — ~ 92, P2) and Vg^l {qi,pi; ~P2, 92) 
in the bound state equation (85), (95) and (109), (117) are determined by the dynami- 
cal mechanism of the effective 4— quark interaction in QCD. In the case of the instanton 
induced 4— quark coupling they are expressed in terms of the form-factors calculated 
in the work of Rapp, Schafer, Shuryak and Velkovskyl 31 ! . For the effective non-local 
4— quark interaction generated by the one-gluon exchange mechanism we have 



Ullt (Pi 



-Qi; -32,f>2j 



4-E 



(pi + qiY 



(pi + qi) 



(7„®Aj)f; (7„®Aj) 



B 2 



a [(pi - pi) 



(pi -P2) 

where a (k) 2 is the running quark-gluon coupling constant in QCD, Xj arc the Gell- 
Mann matrices of the color symmetry group SU (Nc); the flavor indices play no role 
and are omitted.. In reality there do exist the contribution of both above mentioned 
dynamical mechanisms of the effective 4-quark interactions, and they must lie taken 
into account simultaneously. Note that in the special case of the direct 4-quark coupling 
with the interaction Lagrangian (1) or (2) the Bethe-Salpeter equation (5) and (117) 
reduce to the Nambu-Jona-Lasinio equations (38) and (61) resp., This means that the 
Nambu-Jona-Lasinio equation is a particular form of the Bethe-Salpeter equation. 



4 Bethe-Salpeter Equation for Bipolar it ons. 

Now we apply the method presented in two preceding Sections to the study of the 
bound states of two quasiparticlcs in the condensed matters. As a typical example 
of the system of two quasiparticles with the complicated dispersion curves and the 
momentum dependent interaction potential energies we consider the bipolariton-the 
bound states of two excitonic polaritons in semiconductors. The excitonic polaritons 
are the elementary excitations whose formation is the consequence of the quantum 
mutual transition between the photons and the excitons-the exciton-photon mixing. 
Denote 7^ (p) and 7+ (p) the destruction and creation operators of the photon with 
the momentum p and the energy w(p) in the polarization state labeled by the index 
a = 1, 2 , where 

w(p) = ep, 

e being the background dielectric constant of the semiconductor. Due to the elec- 
tromagnetic interaction of the photons with the electrons there arises the quantum 
transition between the photons and the excitons. Applying the second quantization 
formalism to describe the excitons in the local approximation, we denote B a (p) and 
B£ (p) the destruction and creation operators of the exciton with the momentum p 
and the energy E(p) in the spin state labeled also by the index a =1,2 defined such 
that the quantum mutual transition between the photon in the polarization state with 
the index o\ and the exciton in the spin state with the index 02 is allowed if o\ = U2 
and forbidden if o\ ^02- The Hamiltonian of the photon-exciton system without the 
exciton-exciton interaction can be written in the form 

#o = EE{c( P ) 7 ^( P ) 7CT (p) + i?( P )S ff +( P )S ff ( P ) 

a p 

+ ^ [7+ (P) B a (p) + B+ (p) la (p)] J (119) 
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with some function <?(p). The Coulomb interactions between the electrons and the 
holes induce some effective exciton-exciton direct- coupling with the interaction Hamil- 
tonian 



H -°* = EEE< (pi + k ) B £ (pa - k ) u ™ ( k ) B °* (pi) B ^ (Pa) • 

0"1<T2 P1P2 k 

(120) 

The Fourier transform U t7ia2 (k) of the effective interactions potential energies were 
given, for example, in Refsl 22 ' 23 '. The functional integral of the system of interacting 
photons and excitons is 



Z = 



[D 7 ] [D 7 +] [DB] [DB+] 



expiij dtJ2J2 fa M 



i^-«(p) 



7<T (P) + B+ (p,t) 



+ 



ff(p) 



[7+ (P,t) B a (p,t)+B+ (p,t) 7<7 (p,t)] 



(p,t) 
(121) 



ex P 4EEE< (Pi + k ,*) S i (P2 - k, *) U ai „ 2 (k) ( Pl , t) B a2 (p 2 , t) . 



C1CT2 P1P2 k 



If we neglect the direct exciton-exciton coupling, then the functional integral of the 
system becomes 



Zo = J [D7] [Dl + ] [DB] [DB+] 



*^-«(p) 



7ff ( P ) + (p,t) 



+ 



^ [ 7 + (p >t) Bct (p>t) + B + (p,t) 7ff (p,t)] 



S ff (p,t) 
(122) 



In order to describe the bipolaritons we introduce the composite symmetric bi-local 
field 



$<xi<x 2 (Pl,P2,t) = ®<r 2 <n (P2,Pl,t) 



(123) 



depending on two momenta pi and P2 and being labeled by two indices a\, a^- Denote 
$+ i(T2 (pi,P2,t) its hermitian conjugate and set 

Z ** + = J [D$] [L>$+] (124) 
° XP I I / dt EI EI EI *£<x a (Pi + k , P2 k > *) ^kfj (k) *<xi<T2 (Pi. P2, *) > 

K CT1CT2 P1P2 k ) 

Shifting the functional integration variables 

*<Ti* 2 (pi,P2,t) -» *<7 1 a 2 (pi,P2,t) + B ai (pi,t)B a2 (p2,t) ) (125) 
$+ CT2 (Pi,P2,*) -» $+ 1<72 (pi + k,p 2 -k,t) + B+(p 1 + k,t)B+ ( P2 -k,t), 



we establish the Hubbard-Stratonovich transformation 
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exp | - % - f dt e E E B *i (pi + k ' *) B °2 (pa - k < *) ^ ( k ) B ^ (pi> *) B <*> (P2, *) I 

I (Ticr 2 PlP2 k J 

^ I [£>*] (126; 
• exp I - / dt Y E E (P 1 + k ' P 2 k < *) ^ (k) *a lCT2 (Pl,P2,*) > 

I ffl<T2 PlP2 k ) 

.cxpUy^E e t A ^i-2 (pi,P2,t)s CT1 ( P1 ,t)s CT2 (p 2 ,i) 

l CTl<T2 PlP2 



(Pi,i)S+ (p 2 ,t) A fflff2 (pi,p 2 ,t)] }, 



where 



A CT1<J2 (Pi,P2,i) = E^i^ ( k ) *<ria 2 (Pi ~k,p 2 +k,i), 

k 

A+ CT2 (pi, P2 ,t) = E<- 2 (Pl- k 'P2+ k ^)^^(- k )- ( 127 ) 
k 

Using this relation to transform the last exponential in the r.h.s. of the formula 
(121), reversing the order of the functional integrations over the elementary holds 
7 CT , 7+, B ai _B+ and the composite ones $ (Tl<T2 , &t 1 <j 2 and then integrating out over 
the elementary helds, we rewrite Z in the form of a functional integral over the bi-local 
composite field $o- 1<T2 and its conjugate ^tta^ 

Z = T^pr J [£>*] [D$ + ] exp {iS eS [$, $+] } (128) 
with the effective action of the bi-local composite field 

C1IT2 PlP2 k 

$+ ff2 ( P i + k, p 2 k, i) C/ aiff2 (k) $ aiff2 (pi, P2, t) + W [A, A+] , 
where the functional W [A, A + ] is determined by the formula 

exp {iW [A, A+] } = ^ J [Dy] [D 7 +] [DB] [DB+] 



cx P M *EE (p>*) 



<T p 



7a (P) 



( P ,t) 



#a (p,t) 



(130) 



+ ^ [7* (p,*) S CT (p,t) + Bt (p,t) 7ct (P,*)]) } 

.exp^y^E E t A ^- 2 (Pi,P2,t)B«n (Pi,*)S CT2 (p 2 ,t) 
I CT1CT2 P1P2 

+ B+ ( Pl , t) B+ 2 (p 2 , t) A CT1CT2 ( Pl , p 2 , i)] } 



and can be represented in the form of a functional power series in the fields A CTl0 - 2 (pi, p 2 , t) 
and A+ 1(T2 (pi,p 2 ,i) 
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W [A,A+] =Y, w(2n) t A ' A+ ] 



(131) 



n=l 



|y( 2 ™) [A, A + ] being a homogeneous functional of n — th order with respect to each 
type of fields A CTlCr2 (pi,p 2 ,i) and A+ 1(T2 (pi,p 2 ,i) . 

In order to calculate the functional integral in the formula (130) and derive the 
explicit expressions of W^ 2n ^ [A, A + ] we diagonalize the Hamiltonian H of the free 
polaritons by means of the Bogolubov transformation 



B (p) = U (p) Cla + V (p) C 2cr (p) , 
7<r (P) = -V (P) Cla + U (p) C 2a (p) 



(132) 



It can be shown that with the transformation coefficient u (p) and v (p) are determined 
by the relations 



u(p) 

v{p) 2 
2u (p) v (p) 



1 

2< 



1 + 



w(p)-E (p) 



([c(p)-S(p)] 2 + .g( P ) 2 ) 1/2 j 
u(p)-E (p) 



1 

2< 



(Hp)-£(p)] 2 + ff ( P ) 2 ) 

g(p) 

([c(p)-s(p)] 2 + ff ( P ) 2 ) 



(133) 



1/2 = 



the Hamiltonian i?o becomes 



2—1,2 (7 p 



(134) 



where 



Ei (p) 
£ 2 (p) 



i | W (p) + E (p) + (V (P) - ^ (P)] 2 + .9 (P) 2 ) 1/2 } , 

\ ju, (p) + £7 (p) - ([a; (p) - £ (p)] 2 + .g (p) 2 ) 1/2 } . (135) 



Cj CT (p) and cj. (p) are the destruction and creation operators of the polariton with the 
momentum p, energy Ei (p) and polarization a in the branche i. In terms of the mew 
functional integration variables we have 



Z = f [Da] [Dcj] exp ( i f dt E E E °t (P. 



Cta (P, *) 



(136) 



Define the two-point Green functions d (p, t\ — t 2 ) of the free polaritons in the fol- 
lowing manner 



Y Q J [Dei] [Dcf] (Pl^l)C 2 (P2,*2) 

exp|z| ^EEE C ^(P'*) 



l cr p 

iSi 1 i 2 S ai(72 dp lP2 Gi (p^ ti — i 2 ) 



c ilT {p,t)} (137) 
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and denote Gi (p, w) their Fourier transforms with respect to the time variable 



G,(p,i) = ^ / e-^ t G i (p,u)duj 

Z7T 



(138) 



Instead of the time-dependent bi-local fields $,j 1(T2 (pi, p 2 , i), (pi,p2,t) and 

A CTlCr2 (pi,p2,t) , A+ iCT2 (pi,p2,f) it is also convenient to work with their Fourier 
transforms with respect to the time variable $o- 1(T2 (pi,p2,w), &t 1(T2 (pi,P2,w) and 
A CTl<T2 (pi,p 2 ,w), A+ CT2 (pi,p 2 ,o;) 



*<7KT 2 (Pl,P2,£) = ^ j e lUJt ^<r 1 a 2 (Pl,P2,^)dLJ, 

A CTl<T2 (pi,p 2 ,£) = ^ y e" ,ut A fflff2 (pi,p 2 ,w)rfo;. 



(139) 



In the second order with respect to the fields & ai(72 (pi,p 2 ,w), ®t 1 a 2 (Pi>P2j w ) an d 
Acr 1(T2 (pi,p2,w) , A+ 1(T2 (pi,p 2 ,w) the effective action (129) becomes 



SeS $+] 



y^EE A ^ (pi.pa.w) (uo) 

CT1CT2 P1P2 

* CTl a 2 (P1,P2,W) - M(p!,p 2 ,w) A CTlCT2 (pi,p 2 ,w) , 



where 



M(p,q,w) 



2tt/ 



u (p) 2 Gi (p, w - u/) + v (p) 2 G 2 (p, ui - ui') 
u (q) 2 Gy (q, u/) + « (q) 2 G 2 (q, (141) 



Using the formula 



Gi (p,w) 



1 



w - (p) + zO 

and performing the integration over the variables u>, we obtain 



(142) 



M(p,q, W ) = " (P , )2 " (q T, r + U(p)%(q)2 



w - £7i (p) - Si (q) w - Ei (p) - £ 2 (q) 

^(p) 2 ^(q) 2 ^(p)Mq) 2 fl43 s 

w - E 2 (p) - £?i (q) w - £ 2 (p) - -B 2 (q) ' 1 ' 
From the expression (14) of the effective action and the variational principle 

SS cS [$,$+] _ n 



£A+ CT2 (pi,P2,w) 

it follows the field equation 



(144) 



$a 1( r 2 (Pi,P2,w) - M (pi,p 2 ,o;) A CTlCT2 (pi,P2,w) = (145) 
or, in another form 

L (pi, p 2 , w) $ CTl(T2 (pi, p 2 , w) - U aifT2 (k) $ CTlCT2 (pi - k, p 2 + k, w) = 

(146) 

with 
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L(p u p 2 ,u)) = — (147) 

M(pi,p 2 ,o;) 



it becomes the Schrodinger equation for the biexcitons 



[E X (pi) + E 2 (p 2 )] $a lCT2 (Pi, P2, W) + 

+ U*kt 2 ( k ) <JW 2 (Pi k, p 2 + k, w) = w$ fflff2 ( Pl , p 2 , w) (148) 

k 

In terms of the functions of the space-time variables and the differential operators of 
the space-time coordinates we can rewrite the Bethe-Salpeter equation (146) in the 
form 



L (- — - — i— 

\ dr\ ' dr 2 ' dt 



- U air72 (ri - r 2 ) 



$ CTlCT2 (r 1 ,r 2 ,i) = 0, (149) 



where E/<7 1(T2 (ri — r 2 ) are the effective interaction potential energies between two ex- 
citons in the polarization states u\ and u 2 . 



5 Conclusion and Discussions 

In this series of lectures we have presented the universal method based on the applica- 
tion of functional integral technique for the study of the bound states of the systems 
of two relativistic particles in high energy physics or two quasiparticles with arbitrary 
dispersion laws in the condensed matters. The Bethe-Salpeter equation for the compos- 
ite fields describing the corresponding two-body systems were derived. In the special 
case of the direct local 4— fermion coupling the Bethe-Salpeter equation reduces to the 
Nambu-Jona-Lasinio equation. In the derivation of these equations we have observed 
a very interesting phenomenon : there might exist some composite boson fields with 
non-vanishing vacuum expectation values-the dynamical Higgs fields. The existence 
of these fields and the spontaneous breaking of the corresponding symmetries in QCD 
should be studied subsequently. 

On the example of the bipolariton problem we have demonstrated how to apply the 
functional integral method in order to derive the Bethe-Salpeter equation for the bound 
states of two quasiparticles with the complicated energy spectra in the condensed 
matters. The solution of the established bound state equation for the bipolaritons, the 
derivation of the Bethe-Salpeter equations for other two-body systems in the condensed 
matters, for example the biphonons, the bimagnons, the electron-phonon or electron- 
magnon bound states etc.. as well as the solution of these equations would be also, 
the interesting subjects for the subsequent study. 

The presented functional integral method can be easily extended for the applica- 
tion to the study of different three-body systems in high energy physics as well as in 
condensed matters theory. For the definiteness and as an example let us consider the 
formation of the baryons as bound states of the tree-quark systems - the triquark. As 
the necessary tool for the application of the functional integral method with the use 
of the Hubbard-Stratonovich transformation we introduce the six-quark (non-local, in 
general) coupling induced by the interaction mechanism in QCD (instanton induced, 
gluon exchange etc.) with the contribution 



Sint [i>,W[ = \J d A k± J d 4 Pl J Ai J d 4 k 2 J d 4 P2 J d 4 q 2 

(4) (M 

v c?b?m {<H,Pi,k 1 ]k 2 ,P2, q 2 ) ip C2 fe) V>s 2 O2) V% O2) 



S {4) (h +pi+qi-k 2 -p 2 - q 2 ) (h) (px) V> Cl (qi) (150) 
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to the effective action reproducing the 6-point Green functions of the interacting quark 
field. In order to describe the baryons as the triquarks we introduce the tri-local 
composite fermionic field and its conjugate with the Fourier transforms ^>cba {q,P, k) 

^_ABC 

and (k,p,q). We set 

z** = l^^cxpj-^l^ldV/Ai 

J d 4 k 2 1 d 4 p 2 1 d%6 ( ^ (h + Pl + qi -k 2 - P2 - q 2 ) (151) 

— AlBlCl ARC ~ 1 

* VciBiAi (ll,Pl,ki;k 2 ,p 2 ,q 2 )^c 2 B 2 A 2 (<?2,P2,fe) > 

and establish the Hubbard-Stratonovich transformation 

exp jl J d 4 k x J d 4 Pl J d% J d 4 k 2 J d 4 p 2 J d 4 q 2 5 {4) (h + P i+qi-k 2 -p 2 - q 2 ) 

~Ai ~Bi ~d „ 

tp {ki)ip {pi)ip (qi)V Cl 2 Bl 2 Al 2 (q 1 ,p 1 ,k 1 ;k 2 ,p 2 ,q 2 )ip C2 (q 2 ) ip B2 (p 2 ) ^ M (k 2 ) 



= -1= J [DUf] [DW] exp j-1 J d 4 k x J d 4 Pl J d% (152) 
J d 4 k 2 J d 4 p 2 J d 4 q 2 8 {4) (ki + p\ + qi - k 2 - p 2 - q 2 ) 

~A 1 B 1 C 1 ABC ~ 1 

* (ki,Pi,qi)V c 2 B 2 A 2 (q ll p 1 ,k 1 ;k 2 ,p 2 ,q 2 )^c 2 B 2 A 2 (q2,P2,k 2 ) *> 

expj-^ J d 4 k J d 4 p J d 4 q 



~A ~B ~C „ ^_ABC „ „ „ 

tp (k)ip (p)ip (q) Acba {q,P, fc) + A (k,p,q)->Pc(<l)'>pB(p) , >pA(k) 



where 



A CB A(q,P,k) = J d 4 k' J d 4 p' J d 4 q'5 {4) (k+p + q-k' -p' -q') 

Vcb B a C ' (q,p,k;k',p',q')* C 'B>A' (q',p',k'), (153) 
"R BC (k,p,q) = J d 4 k' J d 4 p' J d 4 q'S^ (k+p + q-k' -p' -q') 

f ABC {k',p',q')V£ B B ? A ,{q,p',k'-,k,p,q). 

Substituting the expression (152) into the r.h.s of the formula (4), reversing the order 
of the functional integrations and integrating out over the quark field and its conjugate, 
we rewrite the functional (4) in the new form 

Z=-^= [ [DUf] [Dtt] exp {Seff [*,*]} (154) 

with the effective action 

S eff [*,W] = ~J d 4 k x J d 4 Pl J d 4 Ql J d 4 k 2 J d 4 p 2 J d 4 q 2 
m ^_ A1B1C1 

•<5 W (fci + pi + qi - k 2 - p 2 - q 2 ) # {ki,pi,qi) (155) 

■Vi^M (qi,Puki;k 2 ,p 2 ,q 2 )^c 2 B 2 A 2 (to, to, k 2 ) + W [A, A] , 
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W [A, A] being a functional power series in the composite fields Aqba (q,P,k) and 

^_ABC 

A (k,p,q) ■ In the lowest (second) order approximation we have 

W [A, A] « W^ (2) [A, A] = ~ /" d 4 fc /" d 4 p y A (156) 

A~' 4lSlCl (fc,p, g) (?) (p) < 2 (fc) Ac 2 b 2 a 2 (?,p, fc) . 
In this approximation from the variational principle 

^_ABC J =0 ( 157 ) 

it follows the Bcthc-Salpeter equation for the triquarks 



A CBA {q 7 p,k) = J d A k' J dV / d 4 q'S(k+p + q-k' - p' - q') (158) 

Vcb B a C ' (?,P, *; 9') ^ 2 (<?') SfJ (pO s£ (fc') A C2 b 2 a 2 W,p', k') . 

In the special case of the local direct 6— fermion coupling of the quark field the po- 
tentials Vqba (QjPi k\ k',p', q') do not depend on the momenta k, p, q, k', p', q' and 
the Bethe-Salpeter equation (158) reduces to the Nambu-Jona-Lasinion established 
and discussed in Ref' 32 ' . They would be the basic equations for the theoretical study 
of the baryon structure in QCD. 
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